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ABSTRACT 

We model the X-ray properties of millisecond pulsars (MSPs) by considering hot spot emission 
from a weakly magnetized neutron star (NS) covered by a hydrogen atmosphere. We investigate the 
limitations of using the thermal X-ray pulse profiles of MSPs to constrain the mass-to-radius {M/R) 
ratio of the NS. The accuracy is strongly dependent on the viewing angle and magnetic inclination but 
is ultimately limited only by photon statistics. We demonstrate that valuable information regarding 
NSs can be extracted even from data of fairly limited photon statistics through modeling of archival 
observations of the nearby isolated PSRs J0030-I-0451 and J2124--3358. The X-ray emission from these 
pulsars is consistent with the presence of an atmosphere and a dipolar field configuration. For both 
MSPs, the favorable geometry allows us to place limits on the allowed M/R of NSs. Assuming 1.4 
M0, the stellar radius is constrained to be i? > 9.4 km and R > 7.8 km (68% confidence) for PSRs 
J0030-f0451 and J2124-3358, respectively. We explore the prospects of using future observatories such 
as Constellation- X and XEUS to conduct X-ray timing searches for MSPs not detectable at radio 
wavelengths due to unfavorable viewing geometry. We are also able to place strong constraints on the 
magnetic field evolution model proposed by Ruderman. The pulse profiles indicate that the magnetic 
field of an MSP does not have a tendency to align itself with the spin axis nor migrate towards one 
of the spin poles during the low-mass X-ray binary phase. 

Subject headings: pulsars: general — pulsars: individual (PSR J0030+0451, PSR J2124-3358) — 
stars: neutron — X-rays: stars 



1. INTRODUCTION 

Recent X-ray studies have revealed that a num- 
ber of known rotation-powered millisecond pulsars 
(MS Ps) exhibit predominantly thermal soft X-ray emis- 
sion (Grindlav et al.l l2002l : [ZavlinI 120061 : iBogdanov et all 
|l]06a; Zavlin 2007). The infered effective emission radii 
i?cff indicate that this radiation is localized in regions 
on the neutron star (NS) surface much smaller than the 
expected stellar radius (i?eff «C R) but comparable to 
the classical radius of the pulsar magnetic polar cap 
Rpc = {2TrR/cPy^^R. This finding is in agreement with 
theoretical models of pulsars, which indicate that the 
conditions in the magnetosphere of a typical MSP favor 
heating of the polar caps to ~10^ K by a return flow of en- 
ergetic particles along the open magnetic field lines (sec, 
e.g., Harding & Muslimov 2002; Zhang & Cheng 2003, 
for details). As this heat is restricted to a small fraction 
of the NS, study of the X-ray spectra and pulse profiles 
of MSPs can reveal important information about the star 
such as the radiative properties of the NS surface, mag- 
netic field geometry, and NS compactness {R/Rs, where 
Rs = 2GM/c^). This approach, originally proposed by 
lPa.vlov fc ZavlinI (|1997l ) in the context of radio MSPs, 
can serve as a valuable probe of key NS properties that 
are inaccessible by other observational means (e.g. ra- 
dio pulse timi ng). As shown by P avlov & Zavlin (1997), 
iZavlin fc Pavlovl (|199a) . and [Bogdanov ct al. (2007), a 
model of polar cap thermal emission from an optically- 
thick hydrogen (H) atmosphere provides an excellent de- 
scription of the X-ray pulse profiles of PSR J0437-4715, 
the nearest known MSP. On the other hand, a blackbody 
model is inconsistent with the X-ray timing data and 



can be definitively ruled out. Furthermore, there is com- 
pelling evidence for a magnetic dipole axis offset from the 
NS center (Bogdanov et al. 2007). Finally, the compact- 
ness of PSR J0437-4715 is constrained to be R/Rs > 1.6 
(99.9% confidence). Thus, modeling of X-ray data of 
MSPs appears to be a very promising approach towards 
answering long-standing questions regarding the funda- 
mental properties of MSPs and NSs, in general. 

The prese nt paper represen t s an e xtension of the work 
described in IBogdanov et all ()2007D . Herein we explore 
the detailed properties of the MSP X-ray emission model 
with particular emphasis on its use for constraints on the 
NS equation of state (EOS). The work is organized as 
follows. In §2 we examine the properties of our model; in 
§3 we discuss an application of our model to archival X- 
ray observations of PSRs J0030-h0451 and J2124-3358. 
In §4 we present a discussion and end with conclusions 
in §5. 

2. MODEL ANALYSIS 

The specifics of t he theoretical model em ployed in this 
work are outlined in lBogdanov et all (|2007l ). In brief, the 
model considers a relativistic rotating compact star with 
two identical X-ray emitting hot spots, which in the case 
of MSPs (presumably) correspond to the magnetic polar 
caps. The position of each spot relative to the observer 
is given by the angle tp defined as 

cos ip{t) = sinasinCcos0(t) -f cosacosC (1) 

where a is the pulsar obliquity (i.e. angle between the 
spin and magnetic axes), C is the angle between the pul- 
sar spin axis and the line of sight to the distant ob- 
server, and 4'{t) is the spin phase. We take a non- 
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Fig. 1. — (Left) Representative synthetic hydrogen atmosphere Ughtcurves for a rotating M = 1.4 Mq, R = 10 km NS with two antipodal 
hot spots for the four Ughtcurve classes (I-IV, from top to bottom, respectively), defined by Beloborodov (2002). The dashed lines show 
the individual flux contribution from the two hot spots while the solid line shows the total observed flux. All fluxes are normalized to the 
value corresponding to a face-on hot spot (o = = 0). Two rotational cycles are shown for clarity. {Right) Orthographic map projection 
of the MSP surface for the four pulse profiles. The dashed line shows the magnetic axis while the dotted line shows the line of sight to the 
observer. The hatched region corresponds to the portion of the star not visible to the observer. 

while the total displacement (i.e. impact parameter) of 
the magnetic axis from the stellar center is 

'As' 



rotating Schwarzschild metric as a description of the 
space-time in the vicinity of the NS and include Doppler 
boosting and propagation time delays. This formalism 
is remarkably accurate as long as the spin period of 
the star is > 3 ms where rotation -induced oblateness is 
unimportant (jCadeau et al.l l2007f ) . The NS surface is 
assumed to be covered by an unma gnetized (B < 10^ 
G), optically-thick H atmosphere (jZavlin et al.l Il996t 
iMcClintock et al.ll2004D . The angle-dependent emission 
pattern of this atmosphere differs substantially from the 
isotropic one of a blackbody, which is particularly evi- 
dent in the rotat ion-induced modulatio ns of the X-ray 
flux (see Fig. 1 of lBogdanov et al.ll2007| ). 

Using this model we generate both synthetic spectra 
and pulse profiles (Fig. 1) for the following parameters: 
the effective temperatures and radii of the emission re- 
gion(s) Ti, T2, R\ and i?2, a, C,^ and M/R. Unless noted 
otherwise, we fix the mass to the canonical NS value 
M — lA Mq and vary only R. We also allow for an off- 
center magnetic axis by including offsets in the position 
of the secondary hot spot in latitude and longitude, Aa 
and A(/), respectively, from the antipodal position. The 
corresponding net offset of the secondary hot spot from 
the antipodal position across the NS surface is^ 

As = R cos^"'^ [cos a cos(a-f AQ;)-|-sin a sin(Q;+ Aa) cos A(/)] 

(2) 



Ax = i?sin 



2R 



(3) 



^ Note that equation (12) in lBogdanov et aTl 1 120071) is valid only 
if the angle a is reckoned from the equator towards the spin pole. 
However, as defined by convention, the pulsar obliquity a is mea- 
sured from the spin pole towards the equator. Thus, the correct 
expression for the seondary hot spot offset is given by Equation (2) 
in this present paper. 



We conduct fits to b oth the spectrum and pulse pro- 
file. As pointed out bv lPavlov fc Zavli"iil (|1997D . this is of 
importance because the spectral fits provide tight con- 
straints on Toff and R^s but do not provide useful in- 
formation regarding the system geometry or M/R. Con- 
versely, the pulse profiles provide constraints on the com- 
pactness and geometry of the NS but are less sensitive to 
the parameters that define the emission spectrum (Toff 
and i?off). 

2.1. Extended versus Point-like Emission 

In iBogdanov et al.l (|2007D only point-like hot spots 
were considered, which was sufficient for application to 
the available data for J0437-4715. Here, we also consider 
more realistic extended hot spots in order to ascertain 
the effect of the uncertain physical size and shape of the 
hot polar caps on the observed pulse profiles. To accom- 
plish this, we consider a grid of emission spots across 
the NS surface, which allows us to construct X-ray emis- 
sion regions of arbitrary size and shape. Figure 2 shows 
the resulting pulse profiles for two extreme cases: 1. a 
uniformly heated, circular polar cap of radius 2 km and 
2. a thin annulus (with thickness much smaller than its 
radius) also with a 2 km radius. The value Rpc — 2 
km corresponds to the expected polar cap radius of a 
P ~ 5 ms pulsar. Surprisingly, the physical extent and 
exact shape of the emission regions do not significantly 
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Fig. 2. — Lightcurves for a rotating M = 1.4 Mq, _R = 10 km NS with two antipodal hot spots with values of a and f as in Figure 1. 
The dashed line shows the idealized case of point-like hot spots while the solid line is for the case of a filled circular hot spot {left set of 
plots) and a thin annulus {right set of plots), both of radius 2 km and T^ff = 2 X 10^ K. The bottom panel for each class shows the ratio 
between the two lightcurves shown in the upper panel. All fluxes are normalized to the value corresponding to a = f = 0. Two rotational 
cycles are shown for clarity. 
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Fig. 3. — The location of pulse profile classes I-IV on the a — f 
plane for MSPs with _R = 10 (solid), 12 (dashed), and 14 (dotted) 
km assuming M = 1.4 Mq. See Figure 1 and text for a definition 
of each class. 



affect tlie observed pulse profile, and only differs from 
the point-like case by <1%, except for class IV where 
it differs by as much as ~15%. The largest differences 
occur at the turning points of the lightcurves. 

In reality, due to the rotation of the pulsar the 
polar cap shape is probably di storted (see, e.g., 
iNaravan fc Vivekanandlll983t iBiggs" 19901 with elonga- 
tion or compression in the meridional (north-south) di- 
rection on the NS surface. However, even in such in- 
stances the differences from a circular hot spot with the 
same effective area are relatively small (few %). This 
also means that for emission regions on the NS surface 
comparable in area to that expected for a circular pulsar 
polar cap, the X-ray pulse profiles do not provide useful 
information regarding the details of the region geometry. 
Moreover, the weak dependence of the pulse profile on 
the polar cap size implies that the distance to the pul- 
sar, which is covariant with the emission area through 
the flux normalization (i?^jj/_D^), does not significantly 
affect the pulse profile shape and pulsed fraction. Thus, 
unlike alternative methods for measuring NS compact- 
ness (e.g., using quiescent low- mass X-ray binaries), fits 
to the X-ray pulse profiles of MSPs are not strongly af- 
fected by uncertainties in the distance. Finally, Figure 
2 implies that the pulse profiles are weakly sensitive to 
broad-band calibration uncertainties in the detector ef- 
fective area as well. Narrow-band deviations from the 
true effective area of the instrument (for instance, near 
absorption edges) are also neglegible as the pulse profiles 
cover relatively wide energy intervals (see §2.2). 

Note that the ambiguity in the exact shape of the po- 
lar caps may ultimately limit the accuracy of the con- 
straint on the desires NS parameters using the approach 
described here. Nonetheless, for the observed spectra 
of MSPs (see, e.g.. lZavlinll200Q) . the effective radius of 
the hotter emission region is found to be of order a few 



hundred meters, which compared to the size of the star 
(~10 km) is effectively point-like. Thus, the uncertainty 
in the true shape of the emission region can be overcome 
by considering only the hotter emission component. 

2.2. Limiting Accuracy of M/R Constraints Using 
MSPs 

It has been shown that combined X-ray spectroscopic 
and timing studies of rotation-powered MSPs can poten- 
tially be used to infer M/R of t he NS (iPavlov fc Zavlh] 
Il997t IZavlin fc Paviwl (l998t iBogdanov et all l2007t ). 
Here, we wish to determine whether with further im- 
provement in data quality, M/R can be determined to 
significantly better accuracy. This is essential if this 
method is to be used for reliable measurements of the 
NS compactness and ultimately the NS EOS. In addi- 
tion, it is important to ascertain how the viewing and 
magnetic geometries affect the constraints on the de- 
sired NS parameters. For this purpose, we have car- 
ried out a series of fits to simulated X-ray data of MSPs 
with greatly improved photon statistics than currently 
available. A major advantage of the greatly increased 
data quality is the possibility of high signal-to-noise 
ratio energy-resolved pulse profiles in multiple energy 
bands and phase-resolved spectroscopy for several phase 
windows. As an illustrative example, we consider the 
current version of the proposed effective area curve of 
the Constellation- X observatory X-ray Microcalorime- 
ter Spectrometer (XMS) detector^ and assume a range 
of exposure times, spanning from 250 ks to 1 Ms. We 
generate representative lightcurves for each of the four 
distinct types defined by iBeloborodov (2002) based on 
the visibility of the two antipodal hot spots (see Figs. 1 
and 3)^, assuming both i? = 10 km and M = 1.4 Mq, 
representative of strange equations of state, and i? = 12 
km and M = 1.4 Mq e xpected from c onventional nucle- 
onic equations of state (iLattimer fc P rakash 2001). 

Class I pulse profiles occur in cases when only the pri- 
mary hot spot is observable at all times, while the sec- 
ondary hot spot is in the invisible portion of the star for 
the entire period (see Fig. 1). Note that for two antipo- 
dal polar caps, this can only be achieved for R/Rs > 1.7 
(or i? > 7 km for AI ~ 1.4 Mq) since for a more com- 
pact star the entire stellar surface is always visible. The 
visibility of only a single hot spot causes large uncer- 
tainties regarding the true geometry and especially the 
compactness. This occurs due to the fact that that vis- 
ible hot spot is always found at a small angle with re- 
spect to the line of sight where light bending effects are 
less pronounced, resulting in weak dependence on M/R. 
As a result a given class I pulse profile can be repro- 
duced by a fairly large set of combinations of a, C, and 
M/R. We find that this problem cannot be overcome by 
deeper observations. Indeed, for simulated deep expo- 
sures (~1 Ms) with Constellation-X the 90% confidence 
limits alone encompass the entire range of plausible NS 
radii (7— 16 km for an assumed 1.4 Mq). Therefore, this 
class of pulse profiles are not suitable for tight constraints 
on M/R. 

Pulse profiles that occupy region II in Figure 3 exhibit 

See |http :/ /constellation . gsfc . nasa. gov / 1| 
^ Note that these are defined for the antipodal hot spot case, but 
the regions are similar for offset hot spots for small displacements. 
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Fig. 4. — {Left) Simulated 1-Ms spectroscopic and timing observation of PSR J0437-4715 with Constellation- X XMS assuming R = 12 
km and M = 1.4 Mq. The choice of phase zero is arbitrary. {Right) Best fit confidence intervals for R (assuming M = 1.4 Mq) from pulse 
profile fits to simulated 250-ks and 1-Ms Constellation- X XMS observations of PSR J0437-4715. 




a single broad pulse per rotation. In this configuration, 
the primary hot spot is observable at all times, while 
the secondary is only visible for a portion of the spin 
perio d. The thermal X-ray pulse prof ile of PSR J0437- 
4715 (jZavlin et al.ll2002t lZavlinll2006D appears to be in 
this class. Figure 4 shows simulated observations of PSR 
J0437-4715, assuming i? = 12 km and M = 1.4 M©, 
with Constellation- X XMS and the corresponding best 
fit confidence intervals for R. It is apparent that with 
a substantial improvement in photon statistics, R could, 
in principle, be constrained to better than 10%, and po- 
tentially down to '^2-3% (at 90% confidence) given suffi- 
cient exposure time. Similar constraints are obtained for 
i? = 10 km and M — lA M© as well. Combined with an 
independent mass measurement from radio timing obser- 
vations, this method could lead to unprecendented con- 
straints on the NS EOS. 

Thermal X-ray lightcurves that are found in region III 
(Fig. 2) are characterized by two pulse peaks per rotation 
period. The n earby PSR J0030-f- 451 is a good example 
of th is class (jBecker et al.l l2000t iBecker fc AschenbachI 
I2OO2I see also §3.1). An advantage of these pulse profiles 
is that each of the hot spots provides a dominant con- 
tribution to the flux of one of the two pulses. This fea- 
ture has several practical consequences. First, it may, in 
principle, allow one to determine whether the two polar 
caps are identical in terms of size and temperature. Any 
significant differences would point to deviations from a 
centered dipole, such as displacements of the dipole along 
the axial (magnetic north-south) direction or small-scale 
multipole contributions. Moreover, the "sharpness" of 
the two pulses is strongly affected by M/R (see Fig. 3 



of iBogdanov et aT]|2007t ) resulting in increased sensitiv- 
ity to the stellar compactness. Thus, class III appears 
to be favorable for tight constraints on M/R. We find 
that a 1 Ms simulated Constellation-X observation of 
PSR J0030-I-0451 permits constraints on M/R down to 
a ^5% level (at 90% confidence). 

Finally, class IV profiles are formed when both hot 
spots are observable at all times. As evident in Figure 1, 
these pulse profiles tend to have significantly lower fluxes 
relative to classes I-III since the hot spots are observed 
close to edge-on for the entire spin period. These pulsars 
may not be observable at radio wavelengths due to the 
large impact angle \a — C|, which may exceed the open- 
ing half-angle of the radio emission cone. Such object 
can only be detected in X-rays. At present, identifying 
such objects as MSPs through X-ray timing is difficult 
without a known radio counterpart, although this may be 
possible with future X-ray observatories (see §4.1). These 
factors make class IV MSPs the most difficult to study 
observationally. Most importantly, as seen in Figure 2, 
these lightcurves are much more sensitive to the size and 
shape of the hot spots. This introduces much larger un- 
certainties (roughly an order of magnitude greater) than 
for classes I-III into the spectral and lightcurve fits, re- 
sulting in weaker constraints on M/R. 

Thus, class II and III pulse profile seem to be most fa- 
vorable for tight constraints on M/R. Fortunately, these 
two classes occupy a major portion of the a — C plane'* for 
the plausible range of M/R, implying that niost M SPs 
should fall in these classes. In IBogdanov et all (|2007f ) we 

The same holds true in the cos a — cos ( space. 
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conducted a detailed study of PSR J0437-4715, a class 

11 MSP. In §3, we focus our analysis on two MSPs that 
exhibit class III pulse profiles. 

3. APPLICATION 

As seen in §2, for favorable geometries stringent con- 
straints on MSP parameters are possible. This suggests 
that even with fairly limited photon statistics some useful 
insight into the MSP properties can be obtained. Below, 
we illustrate this through application of our model to X- 
ray observations of two nearby MSPs, PSRs J0030+0451 
and J2124-3358. Along with J0437-4715, these are 
the only thermal MSPs for which the X-ray pulse pro- 
files are of sufficient quality to allow meaningful con- 
straints on NS parameters, especially M / R. The spec- 
tral analysis for each MSP was performed in XSPEC^ 

12 3.0 using the model described above. As shown 
inlPavlov &: ZaylinI (ll99 7D IZavlin fc Pavlov! (|1998D . and 
iBogdanov et al.l ( 2007f ). even the phase-integrated spec- 
tra of MSP are significantly affected by the choice of 
a and C due to the energy-dependent limb-darkening of 
the atmosphere so it is necessary to consider them in 
the spectral fits as well. In order to apply the model 
to the X-ray pulse profiles, it was first convolved with 
the appropriate instrument response, the encircled en- 
ergy fraction was taken into account, and the sky and 
detector background were added. The fit was performed 
by searching the hyperspace for the minimum. We 
fit the pulse profiles by considering the following param- 
eters: the two temperatures and effective radii of each 
hot spot (Ti, T2, and R2) the two angles a and C, 
the stellar radius i?, the offsets of the secondary hot spot 
from the antipodal position (Aa and A0), and the hydro- 
gen column density along the line of sight {Nn). Unless 
noted otherwise, in our analysis we will assume a fixed 
mass of M = 1.4 M© and allow R to vary within the 
rang e of plausible NS radii for this mass (7 — 15 km; see, 
e.g., iLattimer fc Prakashll200l[ ). 

3.1. PSR JOO3O+O45I 

This 4.86-ms solitary pulsar was discovered at ra- 
dio wavelengths in the Arecibo drift scan survey 
(jLommen et al.l[2000l ). It was subsequently detected with 
ROSAT during the final days of this mission. Due to 
the unstable behavior of the failing PSPC detector, this 
data provided no reliable spectral information but re- 
vealed a do uble peaked pulse profile with pulsed frac- 
tion ~50% (iBecker et al.l[2000h . PSR J0030-h0451 was 
revisited in 2001 June 19 b y XMM-Newton for 31 ks 
([Becker fc AschenbachI 12002^. This observation revealed 
a 0.3-2 keV spectrum qualitatively similar to that of 
J0437-4715. In addition, the EPIC-pn X-ray pulse pro- 
file shows a relatively high pulsed fract ion (^50% for 
0.3-2 keV) with two prominent pulses. iLommen et al.l 
(|2000) have estimated from radio polarization measure- 
ments two possible pulsar geometries: a — 8° and P — 1° 
or a = 62° and /3 = 10°, where /3 is the angular separa- 
tion of the line of sight with respect to the magnetic axis 
at closest approach. As the former combination of angles 
cannot produce the observed double-peaked X-ray pro- 
file for all plausible choices of M/R (see Figs. 1 and 3), 
even if we allow for substantial offsets of the secondary 

^ http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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Fig. 5.— {top) XMM-Newton EPIC-MOSl/2 phase integrated 
spectra of J0030+0451 fitted with a two temperature hydrogen 
atmosphere thermal model. The bottom panel shows the best fit 
residuals, (bottom) Chandra ACIS-S and XMM-Newton EPIC- 
MOSl/2 phase integrated spectra of PSR J2124-3358 fitted with 
a two temperature hydrogen atmosphere model. 



hot spots in a and (j), a — 8°, P — 1° can be defini- 
tively ruled out. This result imphes that the observed 
radio interpulse (jLommen et al.|[2000l ) most likely origi- 
nates from the antipodal radio emission cone, not from 
the same emission cone as the main pulse. 

Figure 5 shows the XMM-Newton EPIC-MOSl/2 spec- 
tra of PSR J0030-H0451. They are well described by two 
thermal components. Fitting a hydrogen atmosphere 
model to the spectrum yields Ti = (1.3 — 2.1) x 10^ 
K, i?i = 0.01 - 0.54 km, T2 = (0.3 - 0.7) x 10^ K, 
i?2 = 0.8 - 3.5 km, and Fx ^ 1.5 x IQ-^^ ergs s'^ (0.3-2 
keV), for assumed iVn = (1 - 3) x 10^° cm'^, M = 1.4 
Mq, i? = 7 — 15 km, and all combinations of a and C in 
the range 0° — 90°. The ranges quoted are la limits. As 
expected, the spectr al continuum is fit equ aly well with 
a blackbody model. iBecker fc AschenbachI ([2002) have 
suggested a broken power-law as a possible alternative 
interpretation of th e spect r um of J0030-I-0451. However, 
as pointed out bv IZavlinI (j2007[ ). this model results in 
unrealistic values of Nh and is thus unlikely. 

The fit to the pulse profile (Fig. 6a) was carried out 
by fixing A^ h = 2 x 10^° ^ '^d taking a distance of 

D 300 pc (jLommen et al.ll2006f ). For M = 1.4 Mq, the 
radius is constrained to be i? > 9.4 km or, more gener- 
ally, R/Rs > 2.3 (68% confidence) for all combinations 
of the other free parameters, with a best fit xt — 0.81 
(for 3 degrees of freedom) . Acceptable fits were obtained 
for radii up to ^^20 km. In the case of a blackbody model. 
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Fig. 6.— (a) XMM-Newton EPIC-pn pulse profile of PSR 
J0030-I-0451 for 0.3 - 2 keV. The solid line shows the best fit 
model. (5) XMM-Newton EPIC-pn pulse profile of PSR J2124- 
3358 for 0.3 — 2 keV. The solid line shows the best fit model, while 
the dashed line shows the best fit for a centered dipole (see text 
for best fit parameters). 

wc find that good fits to the X-ray pulse profile require 
implausably large stellar radii (^20 km). As no NS EOS 
mo dels predict radii exceeding ~15 km for AI — f .4 Mq 
(see lLattimer fc Prakash|[2001h . the validity of the black- 
body interpretation is very doubtful. The lower limit 
on the allowed values of a and C is found to be >44'^"''^ 
(68% confidence), obtained when the other angle is equal 
to 90°. Namely if a = 90° then C = 44° or if a = 44 
then zeta — 90°. Finally, the currently available data is 
consistent with both a centered and a displaced dipole 
field. 



3.2. PSR J2124-3358 



nearby [D 

= 4.93 ms (jBailes et al.l 



250 pc) 



PSR J2124-3358 is a 
isolated pulsar with P 

119971) first observed in X -ravs with the ROSAT HRI 



( Becker fc Trlimpeij[T999t) . As the HRI provided no use- 
ful spectral information only a pulse profile was obtained 
with pulsed fraction ^33%. J2124--3358 was later ob- 
served with Chandra ACIS-S for 30.2 ks and with XMM- 
Newton EPIC-M OSl/2 and EPIC-pn for 68 . 9 and 66.8 
ks. respectivelv (IZavlinI 120061 : iHui fc Becked l2006l ). As 
shown bv IZavlinI (|2006[ ). the X-ray emission from J2124- 
3358 is also well described by a two-temperature thermal 
spectrum. Fitting our hydrogen atmosphere model to the 
phase- integrated spectral continuum yields Ti = (1.3 — 
2.4) X 10^ K, Ri = 0.03 - 0.5 km, T2 = (0.3 - 0.8) x 10^ 



K, i?2 = 0.9 - 3.1 km, and Lx = 1.8 x 10^" ergs s'^ 
(0.3-2 keV), for assumed D = 250 pc, TV^ = 1 x 10^° 
cm-2, M = 1.4 Mq, a = 0° - 90°, C = 0° - 90°, and 
i? = 7 — 15 km. The uncertainties quoted represent ±lfT 
ranges. The derived values are consistent with the results 
by Zav lin (2006) . Note th at the diffuse X-ray emission 
detected around this MSP (|Hui fc Beckeill2006f) does not 
contribute appreciably to the point source MSP emission 
as its total luminosity Lx 1 x 10^^ ergs s~^ (0.1-2.4 
keV) is negligible. 

The X-ray pulse profile of this MSP (Fig. 6b) exhibits 
marginal evidence for a faint secondary peak. Given 
that t his feature is evide nt in both the ROSAT PSPC 
(Bec ker fc Trlimpeij [l999t ) and XMM-Newton EPIC-pn 
(|Zavlinll200 6') pulse profile, it is very likely genuine. The 
fit to the pulse profile (lower panel of Fig. 6b) was car- 
ried out by fixing TVr = 1 x 10^° cm'^ (jZavlinl [200^ 
and assuming the dispersion measure derived distance of 
D = 250 pc. Assuming M = 1.4 Mq, the radius is con- 
strained to be i? > 7.8 km (68% confidence), with best 
fit = (for 3 degrees of freedom). The angles a and 
C are constrained to be >12° (68% confidence) for the 
other angle equal to 90° (i.e. if a = 90° then C = 12° 
or if a = 12 then zeta — 90°). Although the suggestive 
asymmetry of the pulse profile hints at the presence of 
an offset dipole, the poor photon statistics permit a cen- 
tered dipole configuration. As with PSR J0030-f0451, 
a blackbody model requires unrealistically large stellar 
radii (>20 km for M=1.4 Mq) for J2124-3358. Thus, 
we conclude that a blackbody model does not provide 
a valid description of the surface emission properties of 
PSR J2124-3358 as weU. 

4. DISCUSSION 

4.1. Searches for Radio Quiet MSPs 

The effect of light bending combined with the (nearly) 
antipodal configuration of the two MSP hot spots ensure 
that the thermal radiation is observable at Earth for any 
combination of a and C,. In contrast, at radio frequencies 
a pulsar is not observable if ja — C| exceeds the open- 
ing half-angle p of the radio emission cone. This brings 
forth the intriguing prospect of detecting and identify- 
ing such radio quiet MSPs in X-rays using blind pulsa- 
tion searches. With the currently abailable X-ray obser- 
vatories {Chandra and XMM-Newton) this endeavor is 
rath er d ifficult due to the intrinsic faintness of MSPs (see, 
e.g.. lCameron et al.ll2007l ) and their relatively low X-ray 
pulsed fractions (<50%). On the other hand, for the next 
generation of large X-ray observatories ( Constellation- X 
and XEUS), the proposed the >10-fold increase in sen- 
sitivity makes such a survey quite feasible. 

Figure 7 shows the X-ray pulsed fraction of a 1.4 Mq 
MSP with a 10, 12, and 14 km radius as a function of a 
and ^, assuming a t hermal emission spectrum like tha t 
of PSR J0437-4715 (lZavhnll2006l : lBogd anov e talEiM) . 
Also shown are lines delineating the region of the a — ( 
plane for which a pulsar with a given radio emission cone 
width is observable at radio frequencies. Note that the 
true opening angle of the radio cone for a given MSP is 
not known and diffi cult to measure relia bly but could be 
as high as -60° (seeE ramer "eFallllQgsl) . If we assume a 
uniform distribution of pulsar obliquities (a) and viewing 
angles (C), for p < 30° a substantial portion (—45%) of 
the MSP population is invisible to us in the radio. On the 
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Fig. 7. — Contours of constant pulsed fraction(soKd lines) for 
X-ray H atmosphere emission from two antipodal hot spots in the 
0.3—2 keV band as a function of the pulsar obliquity (a) and viewing 
angle {Q for a 1.4 Mq MSP with radius 10, 12, and 14 km (from top 
to bottom, respectively). The X-ray emission spectrum is taken to 
be that of PSR J0437-4715. The diagonal dot-dashed, dashed and 
dotted lines show 10°, 30°, and 60° opening half-angles of the radio 
emission cone, respectively. All MSPs found between each pair of 
lines are detectable at radio frequencies (see text for details). 



other hand, if we consider an X-ray timing survey with 
a hmiting pulsed fraction sensitivity of ~10%, only ~5- 
20% (depending on M/R) of the MSPs will go undetected 
as pulsed sources though they will still be detected as 
X-ray sources (provided they arc not heavily absorbed). 
The Galactic population of MSPs may in fact be pref- 
erentially clustered in a certain range of a due to the 
poorly understood effects of the accretion and magnetic 
field reductio n processes during the LMXB phase on the 
NS (see, e.g.. iBhattacharva fc van" den Heuvell [19911 for 
a review). A deep X-ray timing survey of nearby (^1- 
2 kpc) MSPs may, in principle, reveal whether this is 
indeed the case. 

The detectability of the thermal polar cap emission 
from MSPs for all combinations of a and further sug- 
gests that the entire MSP population of certain Galac- 
tic globular clusters (e.g. 47 Tucanae, seejG rindlav et al.1 
l2002tlBogdanov et al.|[2006al : ICameron et al..,2007,) could 
be detected using high angular and temporal resolution 
X-ray imaging and timing (e.g. with the aid of the pro- 
posed Generation-X observatory). Identifying the whole 
population of MSPs in a cluster has important implica- 
tion for studies of g lobular cluster evoluti on and internal 
dynamics (see, e.g.. lCamilo fc Rasioll2005l and references 
therein). 

4.2. Constraining Magnetic Field Evolution Models 

As shown inl Bogdanov et all (|2007| ) and in this present 
paper, the morphology of the thermal X-ray pulse pro- 
files of MSPs offer useful insight into the magnetic field 
geometry of the pulsar (Fig. 1). This has important 
implications for pulsar magnetic field evolution models. 
For instance, Ruderman (1991; see also Chen & Rud- 
erman 1993 and Chen et al. 1998) has postulated the 
presence of crustal "plates" on the NS surface formed 
by shear stresses on the crust caused by neutron super- 
fluid vortex lines pinned to lattice nuclei (see Fig. 3 of 
Chen et al. 1998). The motion of these plates would 
cause the magnetic fields of MSPs to migrate across the 
stellar surface, resulting in either an aligned magnetic 
field or one "pinched" at the spin pole. However, these 
configurations would result in very little (^few percent) 
modulation of the thermal X-ray flux due to the close 
alignment of the polar caps with the spin axis, regard- 
less of the viewing geometry. This is at odds with the 
observed thermal X-ray pulsed fractions of PSRs J0437- 
4715, J0030+0451, and J2124-3358, which are in the 
range 30-50%. Thus, although the model of Ruderman 
(1991) can reproduce the observed radio properties of 
MSPs, it is inconsistent with the observed thermal X-ray 
pulse profiles of MSPs. In particular, the X-ray pulse 
profiles of PSRs J0437-4715, J0030-HG451 and J2124- 
3358 indicate that the magnetic axes of these MSPs are 
significantly misaligned from the spin axis. This, in turn, 
implies that the dipole fields of these (and likely all) 
MSPs do not have a tendency to align with the spin axis 
nor migrate towards one of the spin poles. 

5. CONCLUSIONS 

We have examined the properties of our model of ther- 
mal emission from hot spots on the surface of a neu- 
tron star covered by a hydrogen atmosphere, relevant for 
MSPs. Our investigation has demonstrated that energy- 
resolved modeling of the thermal X-ray pulse profiles and 
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phase-resolved spectroscopy of the continuum emission 
can, in principle, be used to determine M/R and the 
pulsar geometry to high accuracy, given observational 
data of sufficient quality and favorable values of a and 
As shown in §2.1 the thermal pulse profiles are sur- 
prisingly insensitive to the details of the polar cap size 
and shape, the distance to the pulsar, and the uncer- 
tainty in the instrument effective area. This method 
represents the only feasible approach of studying MSP 
magnetic field, surface properties, and compactness, es- 
pecially for isolated MSPs. Barring any deleterious effect 
such as additional hidd en spectral components (see, e.g, 
iBogdanov et al.l [2006bL and references therein) this ap- 
proach can, in principle, lead to unprecedented insight 
into NS properties. 

Our model is found to be in agreement with the 
observed emission from the nearby solitary MSPs 
J0030-H0451, and J212 4-3358. As with PSR J0437~4715 
(jBogdanov et al.ir2007f ). the relatively large pulsed frac- 
tions observed in PSRs J0030-I-0451 and J2124-3358 re- 
quire the existence of a light-element atmosphere on the 
stellar surface and cannot be reproduced by a blackbody 
model for realistic NS radii. For J0030-f 0451 and J2124- 
3358 we are able to place interesting limits on the al- 
lowed compactness of i? > 9.4 km and R > 7.8 km (68% 
confidence) assuming M =1.4 Mq. Based on our find- 
ings in §2, we expect deeper observations of this MSP 



to lead to much tighter constraints on M/R, which in 
turn may firmly rule out certain families of NS EOS. The 
available thermal X-ray pulse profiles also provide useful 
constraints on magnetic field models of MSPs. Specif- 
ically, the positions of the magnetic polar caps on the 
NS surface implied by the X-ray data indicate that the 
magnetic field closely resembles the conventional oblique 
dipole model of pulsars rather than more exotic field con- 
figurations. 

The success of this approach towards elucidating cru- 
cial NS properties motivates further studies of the nearby 
sample of MSPs with both Chandra and XMM-Newton. 
Moreover, this makes MSPs particularly important tar- 
gets for upcoming X-ray mission such as Constellation- X 
and XEUS. The great increase in throughput of these 
facilities will allow searches for new MSPs, detailed ob- 
servations of a larger sample of known radio MSPs, and 
unprecedented constraints on key NS properties, espe- 
cially the NS EOS. 
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